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ABSTRACT 

This thesis presents an interactive naval radar model 
which computes radar detection in the presence of land 
masses, using a parametric terrain description. 
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I* INTRO DUC TION 



This thesis presents a novel naval radar model which 
computes radar detection in the presence of land masses. The 
model is an interactive computer program which accepts 
scenarios and radar parameters from the user and displays a 
map of the area indicating where targets can and cannot be 
detected. The resulting map can be displayed at the user's 
computer terminal or printed offline. 

Major capabilities of the model are: 

1. parametric terrain description 

2. user friendly input and output 

3. beach return masking 

u. radar shadowing 

5. side lobe masking 

The program is written in FORTRAN IV H (extended ) , to be 
executed on an IBM 3033 using an IBM 3278.2 video computer 
terminal. A data file which contains a parametric terrain 
description must be prepared before using this program, but 
all other required input is prompted at the terminal. 
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II . MOT I V^ION 

The Warfare Environment Simulator (WES) used by Command, 
Communications, and Control personnel is a large scale com- 
puterized naval wargame. The land displayed by the system 
does not affect the radar detection probabilities. This is 
one of the artificialities in the game. 

During Exercise BRIGHT HORIZON *81 the tactic of con- 
cealing small vessels in fjords and among islands to prevent 
radar detection until a surprise attack could be made was 
demonstrated to be effective. This technique, while it 
could be anticipated, can not currently be modeled in WES. 

North Atlantic Treaty Organization naval warfare scenar- 
ios virtually always include the proximity of land. United 
States Navy scenarios are primarily open ocean. It seems 
necessary to incorporate NATO problems into USN procedures. 
With this in mind, a method of encoding terrain data was 
combined with a radar model to produce a method of producing 
maps of the area in question with the concealed areas 
displayed. 



10 



III. MODEL 



A. WARFARE ENVIRONMENT SIMULATOR 
1. Basic Svstenj 



WES is a computer program that simulates a naval 
varfare environment in a large computer system, enabling 
personnel to engage in realistic wartime scenarios without 
the expense of actually losing ships and lives [Ref. 1]. 

2. Ra dar E quatio n 

The WES model uses the standard radar equation modi- 
fied to handle ducting and clutter effects. Equation (1) is 
used to calculate signal excess. 



SE=Pt+2G+2W+TCS-4DR-3-NF-L-C (1) 

SE = signal excess, dB 
?t = peak transmitted power, dB//w 
G = antenna gain, dB 
W = wavelength, dB//cm 
TCS= target cross section, d3//m 
D = ducting factor 
R = target range, dB//nm 
3 = receiver noise figure, dB 

NF = noise factor, dB 



1 1 



C = clutter factor, d5 
L = system loss factor, dB 

System loss factor includes antenna pattern loss and 
atmospheric absorption loss. Antenna pattern loss takes 
into account the change in illuminating energy levels 
brought about by the lobe shape of the radar main beam pat- 
tern. The internal system losses are assigned to be 1.5 d3, 
a typical value for mcst radars. Atmospheric absorption loss 
is added to this. Atmospheric absorption loss is frequency 
dependent, and is assigned as 1 dB around 300 .1Hz and 3,5 dB 
at 5000 MBz, at target ranges typical to naval radars. 

The clutter factor describes the losses due zo sea 
cluster. It is taken zo be lOlog(Hw), where Hw is the signi- 
ficant wave height in feet. 

The ducting factor is used to describe the effects 
of ducted propagation, especially surface evaporation ducts. 
An approximate fit of the IRBPS model of a leaky wave guide 
is used.' A strength of zero indicates no ducting, and a 
radar horizon limitation is evoked. If a duct is present, 
rhis restriction is lifted. A strength of five corresponds 



‘Integrated Refractive Effects Prediction System, 
develooed by Naval Ocean Systems Center and implemented on 
the Hewlett-Packard model 9845 desktop calculator. It is a 
shioboard environmental data orocessor and display system 
designed to aid in the assessment of the impact of lower 
atmospheric refraction effects on Naval EM systems. 
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to a "perfect", lossless duct with free space propagation 
within its boundaries. Values for D are selected from Table 
( 1 ) . 



Table 1: Ducting Factor 



stren gth! 0 1 2 3 

D |1TU TT15 TT:2^3 T7TT' 



The noise factor is used to include jamming. It is 
not used in this model. 

See subsection 3 below for the method of processing 
signal excess. The method used in determining the radar tar- 
get cross section of ships is discussed in Chapter III, sec- 
tion c, subsection 3. - 

For this model an additional term is included. The 
"obstacle gain" is calculated and subtracted from the right 
hand side of the equation . 2 



3* Si^n^ Exce ss mode l 

Signal excess is converted to probability of detec- 
tion on the basis of false alarm number and pulse integra- 
tion as discussed below. 



False alarm number is the numner of radar signal 
pulses that are received before a non-signal noise pulse is 



2?or obstacle gain see chapter III, section c, subsec- 
tion 6 . 
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rsceivsd. WES, as well as this acdel, offers the use of four 
different values. 

The number of pulses illuminating the target is 
evaluated by use of equation (2) , 



Sp= 



60 (PRR) (3W) 
(SW) (ARR) 



where : 



(2) 



N? = number of pulses 

PRR = pulse repetition rate, pulses per second 
BW = horizontal beam width, in degrees 
SW = angular width of swepth sector, in degrees 
ARR = scan rate, in scans per minute 

The distribution of signal excess can be closely 
approximated by a normal distribution with a mean defined by 
the number of pulses integrated, probability of false alarm, 
and a standard deviation of 7dB. A large table contains the 
mean of the signal excess distribution (mu) in terms of num- 
ber of pulses integrated and probability of false alarm. 

Conversion of computed signal excess (SE) to prob- 
ability of detection is accomplished in WES as follows. 

The signal excess is normalized by equation (3). 

2 = (se-mu)/sd (3) 

The 2 is used to enter a table of cumulative standardized 
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noraal values to gat the probability of detection which 
is used in the MES simulation. 

This model runs the process in reverse; given the 
desired probability of detection, the number of pulses, and 
an assumed standard deviation the standardized normal prob- 
ability distribution is used to find a corresponding signal 
excess, which can then be compared to the realized signal 
excess to determine if the target will be detected within 
the constraints given. 

Conditions and As sum pr ions 

A number of constants are established in the i#ES 
model, as well as in this modified version. Many of them 
will be discussed in the sections in which they are singu- 
larly appropriate, but some of the more general ones will be 
mentioned here. 

A 4/3 earth radius is used in calculating the dis- 
tance ~o the radar horizon to compensate for standard atmos- 
pheric refraction. 

Receiver noise figure is taken as a constant value 
of 5.5 d3, a value which is primarily thermal noise under 
standard conditions. 

The return from land is not directly modeled. The 
geometrical cross section of a portion of parametrized 
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terrain is assumed to be the radar cross section. Data is 
unavailable on the reflectivity of varying types of surface 
to coapute the true radar cross section given the geometric 
cross section, so a reflectivity of 1 is assumed. A slope 
of at least .02 is required for beach return, when checking 
from the seaward side, and a slope of at least .3 is 
required for return in the sidelobes. In either case, the 
land must be at least as high as the target vessel. These 
values were determined by experimenting with the geometrical 
cross section of a terrain sample. The steeper slope 
required for the sidelobe masking is caused by the lower 
gain present in the sidelobe. 

3. S:iOLSR-hILLS SODSI 
1 . Basic System 

In 1979 Josef Smoler wrote his Master's thesis on an 
"Operational Lanchester-Ty pe Model of Small Unit Land Com- 
bat" [Hef. 2]. It is a time sequenced, deterministic, batta- 
lion-level, f orce-cn-f orce model implemented on a digital 
computer. It contained a method of modeling terrain devel- 
oped by Christopher James Needles in March, 1976 for his 
Master's Thesis "Parameterization of Terrain in Army Combat 
Analysis" [Ref. 3]. In September 1980 Glenn M. Mills 
attempted to overcome shortfalls in the original model, and 
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added several enrichments to provide added user flexibility 
[Sef. 4], A user's manual is provided for this model on a 
permanent disk in the W. B. Church Computer Center, These 
Fortran programs conduct their Army battles over land, using 
Needles' terrain model and Professor Hartman's elevation and 
line-o f-sight routines. ^ These routines were used to model 
the land in this naval model. Only minor modifications were 
required, partially tc adjust to the change in scale and 
partially to remove some strictly land effects. 

Terrain Model 

In 1976 Christopher James Needles presented and 
evaluated a methodology of parameterizing terrain for use in 
land combat analysis [Ref. 3]. This was a shift from the 
tradirional method of digitizing data compiled from terrain 
and interpolating. He described a method by which terrain 
could be created mathematically by using a modified bivari- 
ate normal probability density function. 

The common form of the bivariate normal density 
function is given by equation (4) . 



3See sections three and four of this chapter for ELEV 
and LOS. 
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The equation has been modified by making the normal- 
izing constant equal the the maximum elevation of the 
desired terrain. The resulting equation has sufficient par- 
ameters to model many different types of "hills'', and by 
carefully combining these hills different types of terrain 
can be modeled. 

This method of terrain modeling is currently being 
used by the STA3 model written in SlilSCRIPT II. 5, but the 
Fortran version of the method has also been used, in the 
Saoler-Mills model previously described- ♦ 

See Figure (1) for an illustration on how the hills 
are constructed. 

3* Si eva t ion Routine 

The purpose or the function ELEV is to compute ter- 
rain elevation for given 2, I coordinates. 

The elvation routine was initially provided to Josef 
Smoler by Professor James K. Hartman for use in his land 
combat model, and subsequently by Mills in his modification 
cf Smoler's model [Ref. 5]. 

The point defined by the coordinates is examined to 
determine if any of the "hills" are present, and if so which 

♦Simulation of Tactical Alternative Response, a 
ground-air combat model developed at the Naval Postgraduate 
School . 
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-•"f.shaoe of fhe hills from above is 
moarriea by the correlation coefficient 






ones. Only those which are effective are evaluated, minimiz- 
ing the computer time required for each point. Given the 
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index of the "effective hill", the appropriate paramenters 



are selected and applied to the previously defined modified 
bivariate noraal distribution to derive the terrain eleva- 
tion at that point. 

In a second elevation routine, ELEVG, the gradient 
components of the equation are computed as well as the ele- 
vation, allowing the slope to be determined easily. 

The procedure to add tree height to the top of the 
hills to increase elevation has been removed for this radar 
model. 

U . L ine of Sight Ro uti ne 

The line-of-sight routine was also developed by 
Professor Hartman [Hef. 5]. Its development and use follow 
closely those of the elevation routines. 

Line of sight (LOS) is a purely geometric computa- 
tion assuming perfect visibility. The result of the LOS com- 
putation was initially the percent of the vertical height of 
the target visable to the searcher, and the percent visible 
of the searcher by the target. In addition, the range and 
elevanion of the point which causes the target to be totally 
obstructed is now returned. 

The basic procedure is: "Find the lowest sight line 
from the searcher over the terrain. Extend this line to the 
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target's position and compare its extrapolated height to the 



target's elevation. Thus compute the percent visible." 

In as much as the searcher and the target are both 
ships, the searcher is assumed to be at sea level. If duct- 
ing is present, the target is also assumed to be at sea 
level. If ducting is not present, an additional routine com- 
putes the amounts by which both the obstructing terrain and 
the target are lowered by the earth's curvature. 

As in SLEV, forest calculations have been removed. 

C. AUG.^ENTATION 

"I • In tro du ction 

The STAS LOS and SLEV routines have been combined 
with the VSS radar equation to get a routine to calculate 
detection behind and next to land. 

A separate routine was written to calculate beach 
effects, side lobe effects, and a special routine to pro- 
vide output on an I3fl 3278 computer video terminal. 

After accepting preliminary terrain, radar, target, 
and environmental data the program samples sites a desig- 
nated distance apart from a specified reference point to 
create a matrix containing either elevation or detection 
data,' to be examined at the terminal, with an option to 
print it on a line printer. 
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Figura (2) presents a flew chart of the overall 

t 

program . 

2 • Ro utine s 

The following subroutines were developed to input 
the infornation required by the programs, 
a. Terrain Data 

k number cf different items are required to 
define the terrain. The elevation of the reference plane 
must be entered, along with the total number of '*hills". For 
each hill, the X and Y coordinates, its maximum elevation, 
the standard deviation in both the X and Y direction, and 
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the correla-ion between them must be entered for each hill. 



This defines the hill, in as much as these are the appropri- 
ate parameters in the bivariate normal distribution 
described in the section on the terrain model. ^ 

The data for the terrain is assumed to be con- 
tained in a data file accessable by the virtual machine. Two 
routines adapted from ones writted by Mark L. Yount perform 
the necessary file definitions, after requesting filename, 
filetype, and filemode (filemode defaults to a) at the 
terminal . 

Information for up to 100 hills covering over an 
area 200nm by 200nm may be entered this way. 
b. Radar Data 

Radar information may be entered in either of 
rwo ways: First, it may be stored in a data file, and loaded 
similiar to the method utilized to load terrain data. Sec- 
ond, it may be entered at the terminal, in response to ques- 
tions on the screen. Table (2) lists the information 
required . 



5?or 3VN equation used in the terrain model, see cha? 
ter III, section B, subsection 2. 
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Table 2: Required Radar Data 



NASS 0? RADAR SYSTEM 

(DEGREES) 

S-E COORDINATES 'OF RADAR (NM) 

NUMBER OF LOBES (UP TO 5) 

STRONGEST i LOBES (DEGRESS CW FROM MAIN LOBE) 

(note: rhe ”#'• above represents the number of lobes desired) 

c. Target Data 

The target data is entered in the same subrour- 
ine. Table (3) lists the required information. 



PEAK TRANS 
RECEIVER I 
PULSE REPI 
HORIZONTAL 
VERTICAL 3 
ANGULAR WI 
SCAN RATS 
PULSE LENG 
FREQUENCY 
ANTENNA HE 



MITTED POWER (WATTS) 
F 3ANCWIDTH (MHZ) 
TION EATS (PPS) 

BEAM WIDTH fDEGREES 
SAM WIDTH ^DEGREES) 
DTH OF SWEPTH SECTOR 
(SPM) 

TH (MICR0SECCND3) 
(HERTZ) 

IGHT (FEET) 



Table 3: Target Data 



IDENTIFIER FOR TARGET 
VERTICAL SIZE OF TARGET (FEET) 

TARGET MAXIMUM DISPLACEMENT (KILOTONS) 



d. Other 

Additional information is required for the pro- 
gram. Ir consists of environmental and probabilistic data. 

Table U: Environmental and Probabilistic Data 

DUCTING STRENGTH (INTEGER-0 FOR NO DUCT TO 5 FOR PERFECT DUCT) 
PROBABILITY OF FALSE ALARM 

REQUIRED PROBABILITY OF DETECTION FOE SINGLE SWEEP 
SIGNIFICANT WAVE HEIGHT (FEET) 

Table (4) lists these additional requirements. 



25 



to be 
is an 



A 

printed 
exam pie 



copy of initial conditions is then 
on the line printer upon request, 
of this listing. 



available 
Figure (3) 
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Figure 3: Initial Conditions 




I I I 




Finally, inf or nation relevant to the actual 
plotting is requested. Simply enter the number of nautical 
miles between sample points desired and a reference point on 
rhe overall map and the routine commences its calculations. 

3. Compu ting Tar get Cross Secti on 

The Target Characterization Branch of the Naval 
Research Laboratory has made carefully controlled measure- 
ments of the radar cross sections of a number of naval 
ships, i simple empirical expression has been obtained (see 
equation (5) ) that expresses the cross section as a function 
of the displacement and the frequency. 

tcs=52f^d^ (5) 

Where: 

tcs - Target cross section is square meters 
f = radar frequency in megahertz 
d = displacement in icilotons (full load) 

This expression provides reasonable approximations 
for the microwave band common to most naval radars, and over 
the displacement range of 2000 to 17000 tons [Ref. &]. 

4 . Be ac h Re tur n t as king 

Each segment of land (elevation greater than 0.5 
meters) is checked to determine if it is a beach. This is . 
done in subroutine BEACH by checking every adjacent sample 
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point in the grid and determining if any of them is ocean 
(elevation less than 0.5 meters) . If this is the case, it is 
assumed that a potential target can get next to it. Next 
subroutine 3RTN estimates the geometric cross section of a 
section of beach. This is done by first determining the area 
within one resolution cell. The size of the cell is deter- 
mined by talcing the area between two circles at the given 
range separated by the range resolution and between two 
radial lines separated by the angular resolution. Subroutine 
SL3VG is rhen called to determine the slope of the section 
of beach at the point of concern as seen by the radar. The 
effective area of the resolution cell is determined by tak- 
ing the sine of this slope and multiplying it by the resolu- 
tion cell total area. See Figure (4) for an illustration. 

If rhe effective cross section determined in the 
above manner is greater than that of the target vessel, it 
is assumed that any return from the beach would mask the 
return from the target. 

5. Nain and Side Lobe I nt er fe rence 

LAND is a subroutine which checks the main lobe and 
given major side lobes for the presence of land. A proce- 
dure similiar to BEACH is used to determine the proximity of 
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A is range resolution 
3 rs angular resolution 
C is anaular resolution 




Figure 4: Geometric Area 



land, and once this is established ELSVG is used to check if 
there is sufficient slope for radar return. 

Side lobes are checked by using ELSVG to check for 
land at the same range as the point to be evaluated in the 
direction of a major lobe when the target is centered in the 
side lobe. 

The oresence of land of sufficient height and slope 
in any of the inspected locations indicates that masking 
takes place. 
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S hadov i no 

Af*ar calling the line-of-sigh" routine, the program 
treats the contact in one of three ways: 

a. Case 1 

There is no obstruction. The obstruction gain is 

set to zero. 

b. Case 2 

There is only partial obstruction. The target is 
at the interface between shadow and illumination. The obsta- 
cle gain is set at the log of an eighth, in accordance with 
the basic optical interference technique. [Ref. 7] 

c. Case 3 

The target is totally obstructed, in which case 
the subroutine OBGAIN is called to determine obstacle gain. 

03GAIN calculates the "obstacle gain" created by 
the radio waves being diffracted over the hill tops. A 
"’•cniie edge" obstacle is assumed, and the calculations are 
performed using Picguenard's radio wave propagation equa- 
tions 1^ the hill is not the highest of the 

radar, the hill, and the target it is assumed that the tar- 
get is again at the interface and the same value mentioned 
earlier is returned, otherwise the calculations continue. 
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In any case, the obstacle gain is treated as an 
attenuation and subtracted from the other values in the 
equation used to calculate signal excess. 

• D isp lay S outi ne 

EMA? is a version of THAP modified to allow the use 
of special symbols. THAP has been submitted for the iJONIMSL 
LI3SAEY at SPS. The routine interactively takes an array and 
prepares a contour map for the vicinity of a given point. 

The output is on a terminal or on the printer, at the option 
of the user, interactively taking input from the keyboard. 
Passed arguments: A is the array to be plotted 

N and a are its dimension 



SFN.a is the number of nautical miles 
between points. 

Other arguments are requested at the terminal, the 
program is a modified version of Professor Gilles Cantinas, 
program CTESAP version of OCT 1, 1969, adapted for terminal 
use . 



commend 

printed 
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ately before 
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top 
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the paper printout. 
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Note: On Figure(6), the representing the radar's posxtion is off the screen 



Figure (6) is an example of what appears on the ter- 
The paper listing gives the entire map area, with 
synibology as before. 
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Figure 6: Screen Display 



In the screen display given, the points are one nautical 
mile apart, with the radar located approximately 20nm to the 
northwest of the center of the island. The letters roughly 
centered in the screen correspond to terrain elevation, as 
coded in Figure (2) . The on the left edge correspond to 

a steep beach, where the beach return is stronger than the 
reflection from the target. The ’** to the right define the 
area shadowed by the island. The clear area around it is sea 
level. 
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8 . Rerun with Ch an2JQ 

This routine allows the model to be rerun after 
changing a select number cf the parameters. Figure (7) is 
the screen display which is the menu of allowable changes. 
Simply enter the number corresponding to the desired change 
and then the new value. The routine will offer to print the 
new parameters and then construct another matrix containing 
elevation and detection data for examination. 
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IV. 



g UTUR E 



A. POSSIBLE UTILIZATION OF MODEL 

This nioael has applications in planning and training for 
Naval Operations. 

For training purposes, this model could be incorporated 
into WES or other sea combat models to better reflect detec- 
tion in areas where land is present. For such uses, the 
input and display routines could be trimmed, and only the 
point in question would be analyzed, not an entire map 
array. 

For planning purposes, ship courses could be selected to 
provide maximum radar protection from land shadowing, or the 
inverse, to allow minimum terrain interference for station- 
ing radar guard ships, could be found after parameterizing 
the applicable areas. 

3. FUTURE ENHANCEMENTS 

Parabolic cylinders could be utilized to model the 
refraction effects instead of knife-edges to improve the 
model somewhat [aef. 9]. 

The input and output routines could be modified to take 
better advantage of facilities available (plotting 
terminals, for instance) . 
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A method to find the beach, instead of identifying if a 
particular point is within a given distance of the sea, 
would enable better coastal coverage. 

Radar return is often observed behind obstacles under 
conditions in which simple knife-edge refraction cannot be 
responsible. The causes for these effecns should be found 
and incorporated into the model. 

For use in other models, different subroutines may be 
adapted for the particular model being constructed. Some 
subroutines will have to be entirely rewritten, especially 
those involving beach return which use mapped data already 
calculated, to handle single-point* problems. 
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APPENDIX A 



FEOGFAM LINING 

The following is an entire listing of the radar model. 

It has been specially written to be run on the IBM 3033 on a 
IBM 3273.2 video terminal, the facilities available at the 
Naval Postgraduate School. Appendix B contains and explana- 
tion of each of the subroutines. 
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UBBOUTINE SETUP 
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APPENDIX B 



SUMOUTINE DMECTORY 

The follcwiag subroutines are used in this project: 

A. BEACH 

This subroutine determines if a point of land is adjacent to 
the sea. 

IX and IX are the X and Y coordinates in array A 
A is the map array, dimensioned N by H 
KOAST is the returned value 

3. 3RTN 

This routine calculates the geometric cross section of an 
area of land the size of a resolution cell which may contain 
a target, 

XI, Y1 are the position of the radar 
X2,Y2 are the position of the beach 
RRES is the range resolution 
ARES is ohe angular resolution 
TMAC is the terrain elevation 

ECS is the effective radar cross section of the beach 
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C. ECDRVS 

This subroutine computes the elevation reduction caused by 
the earth's curvature present at a given range. 

SZS is the effective earth radius in meters 
RANGE is the distance traveled along the earth 
DROP is the loss in elevation resulting from the earth's 
curvature 

D. ELSV 

This routine determines the elevation at a given location. 
X,Y are the coordinates of the point 
TMAC is the returned elevation 

E. ELEVG 

This routine is similiar to SLEV, but it also calculates the 
gradient components. 

GX AND GY are the calculated components of the gradient. 

F. GETSE 

This routine processes radar data to compute the signal 
excess returned to the radar from a contact under the speci- 
fied conditions. 

SE is the returned signal excess 
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G. INPOT 

This routine reads in the radar, target, environmental, and 
detection data, either from the terminal or from a data 
file. (the variables are listed in order of appearance in 
the subroutine) 



HD3S 


name of the radar system 


?TH 


peak transmitted power 


3 


receiver bandwidth 


PRR 


pulse repetition rate 


3W 


horizontal beam width 


37 


vertical beam width 


A33S 


angular resolution 


G 


antenna gain 


SW 


angular width of swept sector 


ARR 


scan rate 


PL 


pulse length 


FREQ 


f requency 


HR 


radar antenna altitude 



TMICA radar antenna altitude 
3X,RI coordinates of radar 





wavelength 


HP 


number of pulses illuminating a target per sweep 


TGTS 


identifier for target 
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HT 


vertical size of target 


SIZ23 


vertical size of target 


DISP 


target maximum displacement 


TCS 


target radar cross section 


I5T3 


ducting strength code letter 


IPF 


probablity of failure code letter 


POST 


required probability of detection for a single 




sweep 



SS signal excess required to give a desired PDET 



HW 


significant wave height 


C 


clutter factor 



3L system loss factor 





receiver noise figure 


saiN 


minimum range 


RRES 


range resolution 


Ha AX 


maximum unambigious range 


RHORIZ 


radar horizon 



H. INTHO 

This is the initial pace display. 
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I. K0V2R 

This routine is only used internally by LOS 

J. LAND 

This subroutine determines if a point would be masked by 
rerurn from land in either the main beam adjacent to the 
given point, or by land at the appropriate point in a side 
lobe. 

NASK is the returned value, one if there is masking, zero 
otherwise. 

K. LOS 

This routine calculates the line-of-sigh t in terms of a 
fraction visible for obser ver-t arget pairs. (the variables 
are listed in the order of appearance in the subroutine) 
XA,YA (X3,Y3) are the X,Y coordinates on the field for A 
and 3 . 

T.'IACA (TNAC3) are terrain elevation for A (B) 

T3ICA (TMIC3) elevation of A (E) . 

SIZ2A (SIZ3B) vertical dimension of A (B) . 

LAT03 (L3T0A) indicator variable for LOS calls. 

HHW is the height of the obstructing hill (if any) . 

WRANGS is the range from the radar to the obstructing hill 
(if any) . 
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VISF3A (VISFR3) fraction of SIZEA(SIZEB) which can be seen 



by 3 (A) . 

L. 03GAIN 

This routine calculates the "obstacle gain 
radio waves being detracted over hill tops 



H1 


elevation of the radar 




HO 


elevation of the hill 




H2 


elevation of the target 




D1 


range from the radar to 


the hill 


D2 


range from the hill to 


the target 


RLAMDA 


wavelength 




OBGN 


returned obstacle gain 





M. $PARS 

This routine is by Mark Yount. It splits a 
string down into its separate words. 

N. RANGE 

This routine finds the range between (XI »Y 
a plane. 

RANGE is returned. 



created by the 



double-precision 



) and (X2,Y2) in 
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0 . 



saAP 



This routine prepares contour maps with special symbols to 
be displayed either on the terminal, in response to loca- 
tions keyed in, or on a line printer. EMAP is composed of 
three subroutines: PREMAP, PSTMAP, and PRTHAP. 

A is the array to be plotted, dimensioned N by M 

SFNM is the number of nautical miles between sample 

points on the grid. 

P. SELECT 

This routine defines the data files so that the program may 
access them. 'Written with assistance from Mark L. Yount. 



Q. SETUP 

This routine reads in the terrain data from a predesignated 
data file. (the variables are listed in the order they 
appear in the subroutine) 

NHILLS total number of hills 
BASE surface reference plane 

XC (I) and YC(I) are the x,y coordinates of the center of 
the hill 

PEAK(I) is the maximum elevation of hill I. 

SX(I) and SY (I) are the standard deviations in X and Y of 
hill I. 
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SHO (I) 

LSI 

NHL 

NHTOT 

LISTH 



correlation between X and Y for hill I 
pointer on where to start reading on LISTH 
array holding how many hills are in corresponding 
grid square 

total number of hills in LISTH 
a list of which hills affect which gridsquare 
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